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FISSION- INDUCED PLASMAS 


By 

W. L. Harries 1 and Yuehjaw J. Shiu 2 
1 . SUMMARY 

The purpose of the work was to investigate the possibility of creating 
a plasma from fission fragments, and to utilize the energy of the particles 
to create population inversion that would lead to laser action. Eventually 
it is hoped that the same medium can be used for both fissioning and lasing, 
thus avoiding inefficiences in converting from one form of energy to the 
other. 

The work consisted of both an experimental and theoretical part. The 
experimental part was an investigation of various laser materials which could 
be used for nuclear-pumped lasing. The most likely candidate for a fissioning 
material in the gaseous form is uranium hexafluoride - UFg, and experiments 
were performed to investigate materials that would be compatible with it. 

The work was done at NASA Langley in collaboration with M. D. Williams, who 
had constructed the apparatus, and gave very great help. The results are re- 
ported in section 2. The theoretical part can be divided into three parts, 
which are given in sections 3, 4, and 5. 

One of the central problems in understanding a fission-induced plasma is 
to obtain a model of the electron behavior, and some preliminary calculations 
are presented in section 3. In particular, the rates of various processes are 
discussed. The purpose of this section is to provide a handy compendium of 
processes for reference. Some sections, especially 3.1.1 and 3.1.2 are largely 
due to J. W. Wilson of NASA Langley, but are included for convenience and because 
the numbers are needed later. A simple intuitive model of the electron energy 
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di^trbution function is also shown. The calculations are only intended to 
obtain orders of magnitude in many instances. The results were useful for 
considering a mathematical model of a nuclear- pumped laser. 

Next a theoretical model of a 3 He-Ar nuc lea r- pumped laser is presented. 
The NASA Langley team were the first to achieve volume pumped lasing in a 
nuclear- pumped laser, at the Aberdeen Proving Grounds, Maryland, and hence 
the model could be compared with experiment. The theoretical modeling was * 
done in collaboration with J. W. Wilson, and the results are summarized in a 
paper which was presented at the First International Symposium of Fission 
Induced Plasmas and Direct Nuclear Pumped Lasers, held at the University 
Paris - Sud, Orsay, France on May 23-25, 1978. This paper, "Simplified 
Model of Volumetric Direct Nuc 1 ear- Pumped 3 He-Ar Laser," is incorporated as 
section 4 of this report. The theory showed good qualitative agreement with 
the experimental results. 

A list of talks and paper connected with this grant is given in section 
5. 

The collaboration of NASA/Langley Research Center personnel is 
gratefully acknowledged, namely M. D. Williams and J. Fryer in the 
experimental part, section 2, and J. W. Wilson and N. W. Jalufka in 
section 4. We also wish to thank R. J. DeYoung of Miami University, Oxford, 
Ohio for his help on section 4. 


2. EXPERIMENTAL INVESTIGATION OF NUCLEAR -PUMPED LASER MATERIALS 

2.1. Introduction and Purpose 

Hitherto nuc lear-puraped lasing has been achieved using 3 He (volumetric 
pumping) and lj B and 235 U coatings (refs. 1, 2). The 235 U has the largest 
fissioning energy (165 MeV compared to 0.76 MeV from 3 He) and, in the form 
of UFg, has the potential of achieving a self-critical gaseous lasing system. 
The UFg molecule is not likely to lase, and the gas must be mixed with one 
or more materials which can absorb the fissioning energy to form inverted 
populations. 

However, a number of difficulties could arise. Most of the nuclear- 
pumped lasers are recombination type, requiring a large free electron popula- 
tion in the plasma. The UF 6 is highly electronegative, and could attach the 
electrons to form negative ions, thus reducing the recombination rate and 
quenching the laser output. Moreover, it has many vibrational and rotational 
energy levels and could deexcite many excited species necessary for lasing. 

In this experimental investigation, several laser systems were studied, 
and the effects of UFg were examined in order to find a laser system 
compatible with its presence. The experiments were carried out with 
electrically excited plasmas. 


2.2. Experimental Arrangements 

Two experimental arrangements were used (fig. 2.1). The quartz laser 
tubes were 80-cm long and 2.5-cm outer diameter. The ends were sealed by 
quartz glass windows at the Brewster angle. The electrode arrangements 
were of two types. In figure 2.1a the positive electrode was placed in the 
center of the tube, and electrodes A and B were strapped and grounded. The 
discharge was then along the tube axis or "longitudinal," although the 
fields were reversed in each half. A "transverse" discharge could be 
obtained by placing two brass rods along the axis (fig. 2.1b), but 100 
1-Kft resistors were soldered to rod C, with the wire connectors on their 
other ends facing rod D. Once a discharge occurred from one of the wires, 
a voltage drop occurred along the resistor which limited the discharge 
current, and this has the effect of more even breakdown throughout the tube. 
It also allowed operation at higher pressures. Voltage pulses of up to 
20 kV and of 0.5-us duration were applied to cause breakdown. 



When investigating lasing, mirrors were used as shown (fig. 2.1), with 
the right-hand mirror either a quarter-wave dielectric, or a gold-coated hole 
coupled type. The light output was focused on an indium arsenide diode, which 
was sensitive in the 1.2-u to 3.5-u spectral region, and its signal displayed 
on an oscilloscope. 

The time behavior of the light output is shown in idealized form in 
figure 2.2. Onset of the current usually reduced the applied voltage, which 
explains the "lag" between voltage and current. The laser output during the 
current pulse included the effects of direct electronic excitation. The pulse 
was usually much higher tnan that in the afterglow, and suffered from noise. 
Electronic excitation did not occur in the afterglow laser output, which showed 
a peak about 0.5*us to 1-us after the current ceased and then decayed in a 
few microseconds. 

For spectrographic studies the front mirror was removed and the InAs diode 
was replaced by a variable wavelength filter or by a McPherson monochromator 
with a grating of 300 lines/mm. The output was not entirely independent of 
wavelength and showed a maximum at 3 ym because of the blazing angle of the 
grating. The radiation was then detected by a phototube with an S-l spectral 

O 

response (roughlv 2000 to 10,000 A) or with Polariod film, sensitive from 

o ' o 

''<2500 A to ''<6500 A. 


2.5. Results 

2.3. 1 . CO; . 

(a) C0 2 spectra, virst the spectrum of C0 2 was recorded on 
Polaroid film using the configuration of figure 2.1a. The 
spectrum appeared weak with the exception of two strong 

O O 

bands at 2883 A and 289b A, which were probably due to 
the ionized molecule C0 2 . 

(b) C0 2 - UFg. Next mixtures of C0 2 - UFg (up to 50 percent 
UF 6 ) were tried. The materials appeared to be chemically 
compatible with and without electrical pulses applied. 
However, spectrograms of C0 2 and C0 2 - UFg showed no 
observable difference. No lasing was detected. 
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(c) He - CO 2 . A He-C02 concentration study was done with the 
experimental configuration ~f figure 2.1a, which was elec- 
trically excited at voltages of typically 20 kV. The 
system lased at the 1.45-g line of the 3 l Pj— 3 l Pj transi- 
tion in atomic carbon. A very small amount of CO 2 was 
needed to obtain lasing at 100-torr total pressure. Variation 
of laser output versus CO 2 partial pressure (fig. 2.3) 

shows that the maximum output occurred for a COj concen- 
tration of 0.04 percent. 

(d) He-C02-'JFe. The experiment on He-CC^-UFg mixtures was done 
with fne configuration of figure 2.1b with typical peak 
voltages of 10 kV. Lasing at the 1.45-u carbon line could 
not be obtained in either He-C02 or He-CC^-UFg (up to 10 
percent) mixtures. The causes for not lasing were not 
definitely known; however, the difference in electrode 
configuration and/or possible UF 6 contamination on the walls 
inside the tube may have been the reasons. 

2.3.2. N 2 . - The N 2 -UF 6 spectrographic study was done in both the con- 

figurations of figure 2.1. Type (b) was used to achieve higher gas pressures. 

(a) N 2 spectra. At 5 torr and 150 torr, N 2 band spectra 
appeared very strongly at the second positive transi- 
tion (C 3 ir - B 3 it). Relative intensities are given in 
figure 2.4. The most prominent transitions were the 

O O 

3371.3 A (0,0) and 3576.9 A (0,1). Some of the N 2 
first negative system and CN spectra were also observed. 

(b) N 2 -UF 6 . A total pressure of 50 torr and UFg concentra- 
tions from 1 to 10 percent were used in the N 2 -UFg mixture. 

There were no obvious changes in the emission spectra, 
with the exception of a Si line which appeared at 2881.58 

O 

A with the addition of UFg. No laser action was found at 

O O 

3371 A or 3576 A, possibly because of the slow rise time 
(150 ns) of the electrical excitation (ref. 3). 
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2.3.3. Cl 2 . 

(a) He-Cl 2 . A mixture of He-Cl 2 was found to lase at the 

1 . 587-um Cl line (3d u F9/2-4p “D7/2) (ref. 4). The effect 
of the percentage of Cl 2 on the laser output at a total 
pressure of 200 torr is given in figure 2.Sa. The 
maximum laser output occurred at 0.2 percent Cl?, lasing 
occurred only during the afterglow. 

(b) He-Cl 2 -UFg. Adding UF 6 at a total pressure of 200 torr 
quenched the laser output very rapidly. The laser out- 
put at 1.587 um was decreased by 90 percent with 3 percent 
UFg, and laser action terminated at £4 percent UFg 

(fig. 2.5b). 

2.3.4. Ar . 

(a) He-Ar. The He-Ar (1% Ar, total pressure 100 torr) system 
(ref. 4) lased at the 1.79-um, 2.9-um, 2.2-um, and 1.2-um 
Ar lines. The 1.79-um line (3d[l/2]°i - 4p[3/2] 2 ) was the 
strongest. Laser signals wer* observed during both the 
current pulse and afterglow period. 

(b) He-Ar-UFg. The 1.79-um Ar line was studied with varying 
amounts of UF$ added to the He-Ar mixture (fig. 2.6). 

Curve (a) is the direct electron excitation output, anc. 
curve (b) is for the afterglow. At £ 2 percent UF 6 , 
lasing in the afterglow stopped, while direct electron 
excitation lasing (initial laser pulse) decreased gradually, 
and completely ceased at £ 8 percent UF 6 . 

<..3.5. Xe . 

(a) He-Xe. This system (2®a Xe, total pressure 60 torr) (refs. 
4,5) was found to lase at the Xe lines, 3.5-ym (5d [7/2) ° - 
6p[5/2] 2 ) , 2.02-um (5d[3/2]f - 6p [3/2] 2 ) , 2.65-um (5d [3/2] f - 
6p[l/2)o)* and 3.65-um (7p [ 1 / 2 ] - 7s [ 3/ 2 ] 2 ) • The 3.5-um 
line was the strongest. 
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(b) He-Xe-UFg. The effect of adding UFg to the He-Xe mixture 
lasing at the above-mentioned laser line.* i» shown in 
figure 2.7, which shows lines from the first pulse which 
were excited by direct electron impact excitation. In 
many cases severe arcing occurred due to the large amount 
of UFg in the laser tube, and high percentages of UF6 were 
very difficult to work with. However, laser action was 
observed even up to 20 percent UF$ at 60 to*r of total 
pressure. The afterglow lines are not shown in figure 
2.7, and suffered the greatest attenuation. The 2.65-um 
Xe line (which lased only in the afterglow) was atten- 
uated most, and was completely quenched at < 5 percent 
UF 6 . The 2.02-um and 3.5-um lines were completely atten- 
uated at 7 percent UFg. At a fixed percentage of UFg (10%) 
the laser output vs. percentage of Xe was measured with 
gold-coated hole couu! id mirrors and found to maximize at 
*4 percent Xe (fig. 2.8). 

(c) Ar-Xe. The Ar*Xe system (200-torr total pressure) (ref. 6) 
lased at the xenon lines 2.65 pm (5d[3/2]“ - 6p [1/2 ]q), 

2.62 urn (5d (5/2] * - 6p[5/2] 2 ) . 2.02 um (5d [3/2] ?) - 6p [3/2] j ) 
and 1.73 um (5d[3/2]j - e>p [ 5 / 2 I 2 ) • Of these the 2.65-um line, 
which lased both during the current pulse and afterglow, 
was the strongest, it was observed that all these laser 
lines were from the 5d level of excited Xe. A concentra- 
tion study (fig- 2.9) showed that the afterglow laser 
output had a maximum at about • 4 percent Xe for the 2.65-um 
line. 


2.3.6. Kr . 

(a) He-Kr (ref. 7). A He-Xr mixture (0.1% Kr, 600-torr total 
pressure) lased at the following Kr lines: 

1.784 um (4d[l/2][J - 5p[l/2] 1 ), 2.524 um (4d[l/2]f - 
5p[3/2] 2 ) and 3.0672 um (6p [ 1 . 2] t - 6s [3/2 ] f) . The l.'84-um 
line was the strongest. A comparison of the output for the 
1.78-um Kr line (He-Ar mixture) with the 1.79-um Ar line 


(He-Ar mixture) was made. The total pressure was 600 
torr in each case, and the concentrations were those for 
optimum power, .amely 0.12 percent Kr and 0.6 percent 
argon respectively. Voltage pulses of about 3 kV lasting 
for 0.5 us were applied using the arrangement in figure 
2.1b. The time behavior of the laser pulses was similar, 
with a steep rise to a peak in about 1 us, followed by a 
decay. The peak power of the Ke-Kr laser was about 1.4 
times that of the He-Ar, and the signals dropped to 1/3 
peak value in approximately 1.5 us(Kr) and 0.7 us(Ar), 
with total durations of 2.5 and 1.5 us respectively. 

(b) He-Kr-UFg. With 0.12 percent UFg added to the He-Kr mix- 
ture (0.12% Kr, 600-torr total pressure), the 1.78-um Kr 
laser line, which previously decayed with time in the 
afterglow, was now observed to have a maximum about 1 us 
after the voltage became :ero. The laser duration also 
increased from 2.5 us to 4.5 us. However, higher percentages 
of UFg caused severe arcing, and no definite conclusion 
can be drawn at present. 

2.4. Discussion and Conclusions 

The major purpose of this experimental investigation was to find a laser 
system compatible with UFg for nuclear pumping. So far, most of the laser 
systems under investigation were recombination types: He-Cl 2 , He-Ar, He-Xe, 

etc. It appeared that adding UF 6 to these laser systems resulted in a decrease 
or quenching of the laser output, suggesting that recombination type lasers 
might not be suitable for nuclear pumping. However, among the systems studied, 
the 3.5-um Xe line (He-Xe mixture) exhibited high gain when the excitation 
process involved direct electron excitation from the ground state. Also, the 
He-Xe-UF^ system was observed to lase at £ 10 percent UF 6 with a relatively 
high percentage of Xe ( ? 4°6) for maximum laser output. This system might also 
be promising for nuclear pumping. A mixture of He with 10 percent UFg and 5 
percent Xe at a total pressure of 600 torr is probably the optimum combination 
for testing at the reactor site. 


S 


Not all the above systems however, need necessarily have acted as 
recombination t;-pe lasers, and it is possible that a different pumping 
mechanism could have occurred in the Ar-Xe mixture, which yielded the 2.65- 
um Xe line. In the Hc-Xe system both 7p-*7s and 5d-6p transitions were ob- 
served, whereas in the Ar-Xe system only the 5d-6p transitions occurred. 

The difference in species formation implies pumping of the upper laser leve 

# 

through the formation of Ar 2 , which subsequently transferred its energy to 
the Sd levels of Xe. If so, the pumping would not require electrons for 
recombination, and should work well with UFg. Further studies on He-Xe 
and Ar-Xe should be carried out to provide a basic understanding of the 
pumping process in these two systems. 

A 3 He-Kr nuclear laser pumped by neutrons from a reactor has been 
found to lase at the 2.52-um Kr line (refs. 4,7), but its power output was 
very low compared to 3 He-Ar at the 1.79-um Ar line (ref. 4). In contrast 
here with electrical excitation it was found that the 1.78-um Kr line was 

the strongest of the Kr lines, with the lasing power ratio of * 0.7 

at their peak 1? •. outputs and optimum concentration. Also, the 1.78-um 
Kr line had a ' <ig<;r pulse than the 1.79-um Ar line, so that their total 
energy output was comparable. Our result suggests that the 1.78-um Kr line 
should be further studied with nuclear pumping both by cylindrical and box 
lasers. There should be no optical alignment difficulties because the 1.79 
urn Ar line and 1.784-um Kr line are so close, although mirror adjustments 
seem more sensitive for Kr line than Ar line. 
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(a) Longitudinal discharge. 
Electrodes A and B are 
connected and grounded. 


C 



(b) Transverse discharge. 
Electrode C consisted 
of 100 wires connected 
to 1-K resistors, to 
cause even breakdown 
throughout the vessel. 


Figure 2.1. Experimental arrangements. 
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LASER OUTPUT 
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t (vis) 


Figure 2.2. Idealized time behavior of laser pulses. The leading 
edge of the light signal was subject to noise pickup. 
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Figure 2.4. N 2 second positive spectrum. The relative density was not corrected by the 
quantum efficiency of the phototube. 




He-Xe-UFg 



3.65 urn 




RELATIVE LASER OUTPUT POWER 




RELATIVE USER OUTPUT 



Figure 2.9. Relative laser output power from an Ar-Xe 
mixture vs. Xe concentration (2.65-um Xe 
line) . 
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3. PRELIMINARY CALCULATIONS ON ELECTRON BEHAVIOR IN A FISSION- INDUCED PLASMA 


3.1. Introduction 


The purpose of this section is to obtain a model for the electron behavior 
in a fission-induced plasma and to see which processes dominate. In particular, 
approximate values are obtained for the rates of various processes such as 
recombination, energy loss, ionization, excitation, energy transfer, and 
particle losses. An intuitive treatment of the electron energy distribution 
is also included. We consider the electrons to be in a background of 3 He, 
but the results will also apply for small admixtures of other gases. 

3.2. Production of Electrons 


The electrons are produced .y the products of the action of thermal 
neutrons on He 3 : 

n ♦ He 3 - p(0.S7 MeV) ♦ T(0.19 McV) 

Both proton and triton ionize. On the average, 34 eV are needed to produce one 
electron-ion Dair, and the electron receives an average of 8-eV kinetic 
energy; thus, 20 percent goes in kinetic energy (ref. 1). Hence, the number 
of electrons produced per neutron, Q, is 

Q * [(0.57 ♦ 0. 19)MeV/ (34 ♦ 8)eV] > 2 * 10 4 electrons 

If the fast-burst reactor ptoduces a flux, F, of 10 1 neutrons/cm 2 /s , 
and they are absorbed with an average range r * 6 cm, then the number of 
electrons produced per second per cubic centimeter is 



QFr -1 


3 x 10 20 el/cm ? *s 


20 


A more accurate estimate is based on a calculation by Wilson (ref. 1) 
for the energy deposited/neutron cm 2 , in unit volume, which is 0.3 Me V/ cm 3 • 
neutron * 0 ™- for p « 3 atmospheres. Multiplying this by the flux gives the 
energy input per cubic centimeter per second, namely, 3 * 10 22 eV/cm-*s. Assure 
that at one atmosphere the energy deposited is one-third as much, or 10 22 eV/ 
cm *s. The production rate of electrons is thus 10 22 (eV cra’ 3 s’ 1 ) (54 ♦ 

8) eV * 2 * 10 20 el/cm 3 
assume 

* 2 * 10 2 ° el/cm 3 «s 

' ' prod 

3.3. Loss Mechanisms and Plasma Containment Time 


Experimentally, the pulses are * 1/2-ms long, but the electron contain- 
ment time, t ^ could be much shorter. Rough estimates of can be made 
by (1) assuming the electrons diffuse out of the plasma to the walls, or 
(2) by assuming that recombination dominates the lift time. 

3.3.1. Containment time by diffusion to the walls . - The containment time, 
t , assuming diffusion of electrons to the wall, is given by 


where x is a characteristic distance * 1 cm, and D ( the ambipolar diffusion 
coefficient. The diffusion coefficient for helium has been measured (ref. 2) 
as * 700 t 50 cm*- s" 1 at 1 torr, so * 1 s at 7b0 torr, much longer than 
the observed pulse length. Hence, diffusion is not important. 

3.3.2. Containment time by recombination . - The recombination rate de- 
pends on whether He* or He 2 dominates. For the former, three recombination 

processes can occur, namely (1) radiative, (2) collisional radiative (He* ♦ e 

* ♦ * 

♦ e -*■ He e) , or (3) neutral stabilized (He ♦ e ♦ He ■* He ♦ He). Processes 

(2) and (3) are orders of magnitude faster than (1). However, with the type of 

discharges here it seems He 2 should dominate. In addition to processes 




♦ * 

analogous to (2) and (3), dissociative recombination (He 2 ♦ e •* He ♦ He) 

is also possible. The recombination rates of He* from (2) and (3), and the 

recombination of He 2 , turn out to be the same order of magnitude, which 

yields (dn/dt). for electrons in terms of n, the electron density, 
loss 


3.3.3. Steady-state density of electrons . - Assuming steady state, 
and equating the production and loss rates, yield 

n * 10 1 Vcm 1 

and the containment time, t c> defined by 



is of order 10" s, and much less than the pulse duration of 5 * 10 -4 s. The 
values of n and r are roughly the same irrespective of whether He* or 
He 2 dominates. 


3.4. Electron Energy Loss Rates 

The electrons are created at energies up to several keV, and thereafter 
they lose their energy by inelastic and elastic collisions. 

3.4.1. Energy losses by inelastic collisions . - In the range 20 < e < 

2000 eV, the most rapid energy loss rate for electron^ is due to inelastic 
collisions of ionization and excitation of both neutrals and ions. The rate of 
change of the energy e of one electron with time is 

§£ * -v 1 E l - v 2 E 2 - ... * -vE (1) 

where v is the collision frequency for level n, and E n is the energy 
given up by the electron. The quantity vE is an average for all levels. 


where 


v 

r li 

V ■ T * N ov 
a r 


(2) 


and v r is the electron random velocity, X is the mean free path, N is 
the density of targets, and o the cross section for inelastic collisions. 
Here we are concerned with orders of magnitude only and assume X » constant. 
Then 


v f « (2c/™) 1 / 2 
and 


dt 

dt 



(I r r <' /! 


(3) 


(4) 


Thus the energy of the electrons varies 



where e Q is the initial energy of the 
istic time for inelastic energy loss. 


with time according to 



electron and T. is a character 
in 



e 1 / 2 X /2i 
o 

E 


( 6 ) 


The excitation cross section for He above 19 eV is 0.5 * 10" 7 cirr (ref. 2), 

and the number of atoms/cm 3 at 760 torr is 2.5 * 10 19 . Hence X * (No)" 1 * 

7.4 * 10" 5 m. The value of T. depends on the choice of e , and if 

in r o 

e * 1 keV, and E * 20 eV, then T. * 4 x 10" 10 s. 
o’ in 

In the energy regime from 2 keV down to a few volts, elastic collisions 
will contribute little to energy loss (see the following pages). 
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3.4.2. Energy losses below 20 eV . - Equation (5) will describe the 
energy versus time down to & few eV. The last inelastic collision will take 
place at the first excited state of He, which is 19.6 eV. The impacting 
electron will have an energy somewhat over 20 eV, and will possess * 1 eV 
after the collision. At 19.6 eV the inelastic cross section becomes zero, 
and below this energy the electrons can only lose energy by elastic colli- 
sions. The elastic collisions will be of two kinds--against neutrals 
(assumed at around room temperature) and against other electrons or ions, 
where long-range Coulomb-type collisions take place. 

3.4.3. Elastic collisions with neutrals . - The rate of change of electron 
energy versus time for one electron from collisions with neutrals is 


de 

dt 


-< c v 

e 


(7) 


where k is the average fraction of its energy lost in the collision, and 


k * 2m/M * 1/3000 

for electrons in 3 He. The collision frequency v g here is the elastic 
collision frequency and is equal to v r /* e - From equations (7) and (3) an 
energy versus time relation is obtained when elastic collisions with neutrals 
dominate, 



and if the electron starts with an energy ej then 



-2 

r 1 

>*(f ^ > 

= Cl 

1 * 

V e * / . 


L e,nj 


Here 


T 


e,n 


is a characteristic time to reduce the energy to ej/4: 



- 1/2 

El 


(9) 


( 10 ) 
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Assuming tj * 1 eV and X^, the mean free path for electron-neutral 

collisions in He J at 1 torr is l7.6 * 10“ 3 (ref. 4), then X ■ 2.3 * 10* 7 m 

e 

at one atmosphere, and 


e.n 


io- 


s 


for the range 1 to 1/4 eV. The function e(t) approaches c ■ 0 asymptotical ly . 

The time required to reduce e from cj ■ 1 eV down to 1/40 eV is 24 ns. 

Note T * v" 1 , and the range of e< must be taken into account. 
e,n r * 

3.4.4. Elastic collisions with cold electron:, . - A measure of the rate 
of energy loss of electrons of temperature T to electrons of temperature 
Ti can be obtained from the equipartition time t . Strictly speaking 
this time is a measure of the rate at which equipartition of energy is estab- 
lished between two groups of particles, namely test particles in a Maxwellian 
distribution and at a temperature T, and field particles at a temperature 
Tj. The collisions are assumed to be long-range Coulomb collisions. The 
equipartition time t is defined by the equation (ref. 5) 


dT 

dt 



(ID 


and 


T 

eq 


3.3 « 10 s e*/ 2 
n In A 


(s); Ti -*• 0 


(12) 


with e * T (eV), electron density n in cm* 3 , and A is the ratio of the 
Debye length to the impact parameter. Assume here that T; * 1/40 eV, ;< T, 
and let T be 10 eV and 1 eV respectively. The values of In A change little 
with n and T, and the values for t are shown in the following table. 
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Values of t 


in s for 


Ne cm* 3 In A T » 1 eV T ■ 10 eV 


10 u 

10.7 

3.0 x 10* 7 

m 

1 

O 

1012 

9.6 

3.4 x 10* 8 

10* 6 

10 l 3 

8.5 

3.8 x 10" 9 

10- 7 

10 14 

7.3 

4.4 x 10- 1 0 

1.4 x 


For n ■ 10 13 /cm 3 , T[ > 1 eV, the time for the hot electrons to be cooled 

by electrons at room temperature is of order T g n - For energies less than 

1 eV the electron-electron collisions dominate as x « v 3 . 

eq r 

3.4.5. Elastic collisions with ions . - The ions, being of equal mass with 
the neutrals, should be quickly cooled, and should then cool the electrons. 
However, the equipartition time for the electrons to be cooled by ions at 
room temperature, namely x^ ., is far greater than the electron-electron 
cooling time x (ref. 5): 


T eq,9-e T eq,e-i 


m/M. = 1/6000 


(13) 


where NL is the mass of the ion, and the cooling of electrons by ions can be 
neglected . 


3.S. Electron Energy Distribution Function 

3.5.1. Initial electron energy distribution .* - The fission fragments 
ionize the He 3 and release electrons. For high-energy fragments whose velocity 
exceeds that of the orbital velocity of the electron, the binding energy of 
the electron can be neglected, and the collision regarded as Rutherford 
scattering. The differential cross section subtended is proportional to E^" 1 , 
where E^ is the energy of the fission fragment, and the electr-ns are ex- 
pected to be in an initial energy distribution function, f j ( e ) « e* 2 , where 
e is the energy of the electron. The maximum amount of energy handed over 


This description is based on discussions with J. Wilson of NASA/Langley 
Research Center. 
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to the electron is given by e c * (2m/M)E^, where m is the electron mass, 

and M the mass of the fission fragment. The initial distribution function 

for electrons created by a fragment of energy E f is shown in figure 3.1 

* # * 

(solid curve). The picture fails around 1.5 e^, where is the ionization 

energy, or around 30 eV, as there are no more excitation levels, and the 
curve goes to a maximum and falls off at lower energies. 

As the fission fragment slows down to an energy E f , the cutoff drops 

• 9 • 

to e , and the distribution is now as shown by the dotted curve. The 
C , 2 

distribution function for all the electrons produced by one fission fragment 
is the sum of such curves, and will have a slope steeper than e" 2 . 

3,5.2. Overall description of the average energy of one el pc cron versus 

time . - Consider one electron created at a kinetic energy * 1 keV. Its energy 

first decays rapidly by inelastic collisions with a characteristic time T^ , 

and by elastic collisions with a characteristic time T (T. , T « v" 1 ). 

e,n in e,n r 

In addition, cold electrons can cool the hot ones. Strong electron-electron 
cooling does not take place until the energies are low, as x « v 3 . However, 
at e « 1 eV, t ** T , although there are 10 6 times more neutrals than 

electrons. The characteristic times for cooling are much less than the plasma 
containment time. 

The electrons can gain energy from other electrons when at low energies 
by Coulomb collisions (section 3.2.4), and by collisions of t^e second kind 
with excited neutrals (superelastic collisions). However, the density of 
excited neutrals is probably of the order of the electron density and much less 
than that of unexcited neutrals, so collisions of the second kind are 
neglected . 

The electrons, then, start off in a distribution as in figure 3.1, and then 

continuously lose energy until at low energies. For energies below 1 eV and 

n =* 10 13 , x << t , the containment time, and therefore at low energies 

eq,e-e c 

a Maxwellian distribution of electrons is attained, with a temperature approx- 
imately that of the neutral gas. Note that only electron-electron collisions 

(t ) lead to a Maxwel’ian distribution, while electron-neutral elastic 

v eq, e-e' 

collisions (T ) need not. The latter reduce the electron energy; the former 
e , n 

redistribute it. This model enables an estimate of the steady-state electron 
energy distribution to be made. 
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3.6. Simple Intuitive Description of the Steady-State Electron 
Energy Distribution Function 

The intent of this section is to provide simple intuitive arguments to 
indicate the essential features of the electron energy distribution function 
and compare them with that of Hassan (ref. 6). The duration of a pulse is 
5 * 10"* 4 s, and the lifetime of an electron » 10" 7 s; hence we assume there is 
a continuous constant source of electrons for the duration of the pulse. 

Each electron then starts at a high energy, loses energy with time, becomes 
Maxwellian at low energies, and is recombined. 

The density, dn, of electrons in the energy interval e to e ♦ de 
is proportional to the amount of time, At, each electron spends in this 
interval , 

dn « At 

Now for one electron, 

-1 

limit At _ /de\ 

Ae-K) Ae \dt/ 


where is the rate of change of energy versus time for an average electron. 

Therefore, 

-1 

dn = f(e)de <* de 

where f(e) is the steady-state electron energy distribution function, and 



We now construct f(e) versn e on a logarithmic plot (fig. 3.2). Starting 
from high energy, the region DE shows a sharp drop as e increases, be- 
cause the initial distribution of electrons is the sum of curves similar to 
figure 3.1; the exact shape of DE is not discussed further. Below about 
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35 eV, the electron energy loss is mainly from inelastic collisions, where 
dc/dt « e 1 / 2 from equation (4); hence in CD, f(e) « e"*/ 2 . The lower 
limit of the region C is around a volt (see section 2.2.2). Near C there 
is less and less contribution from excitation, and in BC electron-neutral 
collisions dominate. Here equation (8) indicates de/dt « e 3 / 2 , so f(e) « 
e -3 / 2 . Below B electron-electron collisions dominate, and as t « t , 
a Maxwellian distribution is formed, with its peak at the temperature of the 
neutral gas at some hundreds of degrees Kelvin. 

The values of the ordinate in figure 3.2 are not absolute, and an order 
of magnitude calculation enables a comparison to be made of the average 
density, n A0 , of electrons per unit energy in AB versus n CQ in CD. 

Assume AB extends over a range Ae AR * 1 eV, and CD over a range of Ae^p * 
30 eV. Then, 


\b a 

t AB AB 


CD 

C CD 


Ac 


CD 


where the times t^g, t CQ correspond to those an average electron spends 

in AB and CD. Now t. * t * 10 -7 s, the recombination time, and t~ * 

AB C LU 

T. * 0.4 x 10 -10 s. Therefore 
in 


^AB 

n CD 


10 3 or 10 4 


The relative scale on the ordinate of figure 3.2 conforms with this ratio. 

The essential details of figure 3.2 agree with an electron energy distri- 
bution curve of Hassan (ref. 6) derived from the Boltimann equation (figure 
3.3). Hassan actually plotted ffej/e 1 / 2 , and the curve for f(e) is due 
to Meador and Weaver (ref. 7). The points B and D in figure 3.2 occur 
in the same energy range as figure 3.3. Between B and D Hassan' s curve 
shows f(e) 5 e” 1 , whereas figure 3.2 shows f(e) <* e -1 / 2 near D, and * e“ 3 ^ 2 
near B. Both curves show a Maxwellian distribution to the left of B with 
about 10,000 times more electrons in the Maxwellian distribution below 1 eV 
than there are in the non-Maxwellian distribution above it. 
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Figure 3.1. Initial distribution function fj(e) for 
electrons; the curves are not normalized 







4. SIMPLIFIED MODEL OF A VOLUMETRIC DIRECT NUCLEAR- PUMPED 3 He-Ar LASER* 


4.1. Introduction 

Recently successful lasing was accomplished by the NASA/Langley Research 
Center team in a volumetric nuclear-pumped 3 He-Ar laser (refs. 1,2). It is 
the purpose of this paper to develop a simplified model for the laser, to 
establish which processes dominate and, in particular, to clarify the 
physics of the reactions. 

The physical model (ref. 3) is described in section 2.2, and the source 
terms arising from the 3 He(n,p) 3 H reaction as well as the rate equations 
are given. Steady-state solutions are derived for the densities of the 
species versus the fraction of argon. The dominant pumping mechanism is 
collisional radiative recombination of the argon atomic ions, and subsequent 
radiative cascading into the upper laser level. The reactions following 
recombination are described in section 4.3, and an expression is derived for 
the difference in densities of the argon 3d and 4p levels (a quantity that 
should be nearly proportional to the output of the laser). Comparison with 
experiment indicates good agreement. 

4.2. Physical Model 

Neutrons from a reactor are thermalized and enter a vessel containing a 
3 He-Ar mixture at a total pressure of order 1 atmosphere. The laser excitation 
energy is derived from the 3 He(n,p) 3 H reaction. The proton and triton 
possess 0.57 MeV and 0.19 MeV respectively, and excite and ionize the gases. 

The major fraction of the energy released produces helium ions, He*, 
and metastable atoms. He*. The He* is rapidly converted to He 2 , which in turn 
is lost by charge transfer to neutral argon to form Ar* . The He* is par- 
tially converted to He£, and both He* and He£ ionize the argon by Penning 

* Prepared by W. L. Harries of ODU and J. Wilson of NASA/LaRC and pre- 
sented at the First International Symposium of Fission Induced Plasmas 
and Direct Nuclear Pumped Lasers, held at the University Paris - Sud, 

Orsay, France, May 23-25, 1978. 
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ionization. Collisional-radiative recombination and argon molecular ion 
formation reduce the Ar* population. On recombination to neutral argon, the 
cascading excited atoms can form a population inversion between the 3d and 
4p levels, which leads to lasing. It is of interest therefore to estimate 
the Ar* population density, especially as it relates to the recombination 
rate. The laser output depends on the Ar* recombination rate and various 
branching ratios for the processes occurring after recombination. 

In estimating the Ar* density, rate equations for the various species 
were obtained (ref. 3). The species considered were He*, He;> (the He* 
converts rapidly to He-i, and both are lumped together as N*), He*, and He£ 
(lumped together as N*), Ar*, Ar 2 , and the electron density; thus there 
are five coupled equations. The various reactions are given in tables 4.1 
to 4.3, section 4.5, together with the values of the rate coefficients ob- 
tained from the literature. The rate coefficients are also labeled on figure 
4.1. The densities of the various species depend on the relative argon 
concentration, which is described by the fraction c, where 0 < c < 1. 

The rate equations, with the various processes, labeled in figure 4.1 

are 


^ r - «= S* (He) ♦ xN* 2 - (s ♦ uN Q )N*n - tn 2 N* - (q ♦ rN Q )N*cN o (1) 

^ » S* (He) ♦ (s ♦ uN Q )N*n ♦ tN*n 2 - 2xN* 2 - yN*n - pcN Q N* (2) 


4^- ■ S* (Ar) ♦ pcN N* ♦ (q ♦ rN )N*cN - uN Ar*n - t'Ar*n 2 
dt r o o o o 


- [v(l - c) ♦ wc] N 2 cAr* 

(3) 

dAr 2 , «. + 
j ■ = [v(l - c) ♦ wc]N 2 cAr - S 2 Aron 

(4) 

n * N* ♦ Ar* ♦ Ar 2 

(5) 
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Here the quantities naming the species correspond to the densities, and 
n is the electron density. Where more than one symbol is assigned a 
reaction in figure 4.1, the different rate coefficients correspond to 
different third bodies and are described later. 

The source terms in the equations, S terms corresponding to the rate 
of production of excited and ionised species, are proportional to the rate 
of production of protons and tritons, and therefore proportional to F, 
the flux density of thermalized neutrons. The number of nuclear reactions 
per unit volume depends on the total gas pressure p and on d the tube 
diameter, since the neutron flux generally undergoes appreciable attenua- 
tion in crossing. With p ■ 1 atm, d * 2 cm, about 25 percent of the neutrons 
are absorbed in the gas (ref. 3), and therefore the source terms here are 
approximately proportional to Fp. The terms also depend on c, the Ar 
fraction, and with the above conditions are (ref. 3) 

(He) * 2 * 10 20 (1 - c) 2 /(l ♦ 3c) (6) 

S* (He) * 0.8 x io 20 (1 - c) 2 / (1 ♦ 5c) (7) 

S*(Ar) * 1.2 x 10 21 c(l - c)/(l ♦ 3c) (8) 

Although computer graphical solutions of the densities were readily obtainable 
for the steady state, they did not clearly indicate which processes dominated. 
Therefore, approximate analytical solutions were attempted. The method was 
to use the computer solutions for comparison and to neglect small terms. The 
computer sclutions showed that n, the electron density, was almost independent 
of c as long as c did not approach unity. It might be expected that this 
would be roughly so, because the energy level schemes of He and Ar are fairly 
similar and the recombination times of the different ion species are approxi- 
mately the same at equal ion densities. Thus a trial value of n independent 
of c was assumed 

The magnitudes of the various terms were compared for a neutral density 
» 3 » 10 19 cm’ 3 and (from exact computer solutions) n * 10 13 cm" 3 and 
4 x 10 12 < N* < 10 14 . It has been found experimentally that lasing ceases if 
c is more than 0.2, so c was assumed between 0 and 0.1 when comparing 
terms . 
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The recombination terms s, u, t, can be neglected in equation (2), 
as well as superelastic collisions, y. In the steady state the equation 
then reduces to a quadratic, but it is convenient to give solutions in two 
regimes depending on c: 


[S*(He)/2x] 1 /2 • 

o 

A 

o 

• 

(9a) 

S*(He)/pcN Q ; 

c > 10" 4 

(9b) 


The dominant source term S* corresponds to electrons exciting neutral 
helium atoms. When there is little argon present, the source term is balanced 
by losses of He* to He 2 [eq. ( q a)]. When there is appreciable argon present 
the losses of He* and He* are by Penning ionization of argon [eq. (9b)]. The 
quantity pc.N Q has the dimensions time" 1 , so we define a characteristic 
time for the loss of He* and He? by Penning ionization: 

t* * (pcN ) _1 (10) 

Hence equation (9b), which applies to the range of experimental interest, 
can also be expressed as: 

N* « S*(He)t* ; c > 10" 4 (11) 

The formation of He 2 from He* ij a rapid process, and the quantity N*, 
the sum of the two densities, is essentially the density of He^. The solution 
of N* from equation (1) depends on whether c is less than or greater than 
10" 4 . If c < 10" 4 , the charge transfer terms (q,r) can be neglected. The 
dissociative recombination term s, and the recombination with He as a third 
body, u, are also small. Assuming n = N*, then as N* is known, equation (1) 
becomes a cubic in N*. 

On the other hand, if c >> 10" 4 , the charge transfer terms are now im- 
portant. Again the rtcombination of He^ to He* (s, t, u) can be neglected. 

The solutions in the different regimes of c are 
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p*(He) <• S*(He)/2l 


1/3 


< 10 -4 


(12a) 


N* » < 


S»(H.) ♦ x(s-(H«)/pcN o ) S » (H .) . 

(q ♦ rN^JcN^ (q ♦ rN JcN ’ C > (12b) 


o' o 


o o 


The N + is formed by direct electron ionization of He and by the forma- 
tion of He 2 from the He* ♦ He* reaction. When there is little argon present, 
the He 2 is lost by recombination, (t term), where the dominating process is 
He 2 ♦ 2e -+ He* ♦ He ♦ e. With appreciable argon present, the formation of 
He 2 from He* ♦ He* is unimportant as the He* is depleted by Penning ionization 
of argon. The loss of He 2 is now mainly by charge transfer to argon (q and 
r). 


The quantity (q ♦ rN Q )cN Q has the dimensions of time"l, so equation (12b) 
can be rewritten 


N* * S* (He) t + ; c > lO -4 


(12c) 


where t* = ( (q ♦ r!\ JcNq]” 1 is a characteristic loss time for He 2 by charge 
transfer. 

Equation (3) gives two solutions for Ar + depending on whether c >> 10* 4 
or not. It is convenient to consider large c first. None of the terms can 
now be neglected, and 


Ar 


[v(l 


S* (He) f S* (He) ♦ S*(Ar) 

- c) ♦ wc]N^c ♦ [u'N Q n + t'n 2 ] 


c » 10 -4 


(13) 


where the terms in the numerator are obtained from equations (9b) and (12b). 
The production of Ar* is due to direct ionization through S*(Ar), to Penning 
ionization which depends on S*(He), and also to charge transfer, which depends 
on S*(He). The sum of the three terms can be defined as a total source term 
S. The denominator represents the loss processes: the first parenthesis 
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corresponds to the loss of Ar* into Ar: (v, w terms), the secmd due to Ar* 
recombining (u', t'). Characteristic loss times for these processes may 

be defined: 


t 2 • 

[(v(l - c) ♦ wc^ N*cJ** 

(14) 

*3 " 

(u'N n ♦ t 'n* |" 1 
1 o 1 

(IS) 

and if both 

processes occur simultaneously, the overall loss time tj 

is 

given by 



*r l 

- t 2 - 1 ♦ t 3 -» 

(16) 

Equation (13) can also be expressed as 


Ar* « 

» Stj ; c » io- 4 

(17) 


If c < 10* 4 , the v and w terms corresponding to the formation of 
Ar* can be neglected, and from equation (3) 

S*(Ar) ♦ pcN N* ♦ (q ♦ rN )N*cN 

Ar * ; c < 10 " (18) 

u ' N n ♦ f ' n 
o 

The numerator is the sum of three source terms as before. However, the values 
of N* and N* are those corresponding to c < 10 -4 or given by equations (9a) 
and (12a). The three source terms may be added to give a total source term 
S^, corresponding to small c. In addition, now n ■* N*. The denominator 
represents loss of Ar* by recombination, so 

Ar* » S t, ; c < 10- 4 (19) 

o 

Over the range of c: 


S T 3 
0 3 

; c < 10- 4 

(20a) 

Stj 

; c » 10- 4 

(20b) 
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The density of Ar* is now immediately obtained from equation (4). When c is 
small the Ar 2 density is much smaller than for the other species. Otherwise: 

; c » 10- 4 (21) 


The s 2 n corresponds to loss of Ar 2 by recombination to Ar*, and into the 
4p level, so if a characteristic loss time is defined 

t 4 * (s 2 n ) _1 (22) 

then equation ( 21 ) becomes 

Ar 2 ’ S, “( t 2 ♦ I 3 ) (23) 

The quantity T 3 /(t 2 ♦ T 3 ) is the branching ratio, or the fraction of Ar* 
that recombines through the Ar' channel in figure 4.1. A comparison of the 
simplified expressions for the densities versus c (dotted) with exact 
computer solutions including all the terms (fig. 4.2) shows good agreement 
within a factor of two. 

4.3. Conditions After Recombination 

Laser action results when the upper neutral Argon 3d level population 
exceeds the 4p level. After recombination, the 3d level is fed by the Ar 
neutrals cascading down in energy from the Saha region, Ar' in figure 4.3. 

Some of the atoms may not reach the 3d level because of collisional deexcitation, 
and/or spontaneous emission of quanta to energies below the 3d level. Also, 
the energy levels in the Ar' region exceed those for the formation of Ar 2 , and 
there is the possibility of the reaction Ar' ♦ Ar -*■ Ar 2 ♦ e, the Hornbeck- 
Molnar effect. We have included the effect of collisional deexcitation, 
spontaneous emission, and the Hombeck-Molnar effect by describing them with 


Ar 2 


s 2 n 


1 ♦ 


u'N n ♦ t 'n 2 
o 

(v(l - c) ♦ wc^N 2 cJ 
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characteristic times and including them in our calculations. The resulting 
gain does not agree with the observed laser output if these effects are 
significant, we believe that these effects (shown dotted in fig. 4.3) are 
unimportant. 

The Ar* can also be lost by formation of Ar 2 (fig. 4.3) with the 
characteristic time 12 - The Ar 2 can dissociatively recombine to give a 
neutral argon in the 4p level, which should reduce the laser output (ref. 4). 
We have investigated the possibility that a fraction of the Ar 2 density 
fills levels other than the 4p, i.e., either the 3d level or levels below 
the 4p, and again find disagreement with experiment. Thus the flow paths 
that dominate are those shown by solid lines. Then the difference in 
densities of the 3d and 4p levels can be estimated in terms of the times 
t 2 , t 3 and the lifetimes t,., x of the levels.. The difference in 
densities should be proportional to the gain G of the laser, and 


G ‘ Ar> [ T 3d T3 ‘ l " T up (T2 ‘ 1 * T 3" l) J 

(24) 

r / t 2 \ i 


« S T,. (———)- T I 

(25) 

3d \x 2 ♦ x 3 / upl 


4.4. Comparison with Experiment 


Figure 4.4 shows G and Ar + plotted versus c and compared with the 
experimental results (refs. 1, 2). Error bars are not shown, but the 
experimental results are within ±30 percent. The ordinate values are not 
absolute, and G was normalized to the experimental results. The experimental 
cutoff at c - 0.2 due to dissociative recombination of Art into the 4p level 
occurs when there is still ample density of Ar*, and seems consistent with the 
value of G. 

In summary, for low values of c the c 'tput increases with the fraction 
of argon, but if c is made too high, laser action is curtailed by the forma- 
tion of Ar 2 which dissociatively recombines to fill the 4p level. 
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4.5. Helium/Argon Rate Coefficients 




Table 4.1. 

Helium 

rate coefficients. 



Reaction 

Symbol 

Rate Coefficient 

Reference 

HeJ 

♦ 

e He* ♦ He 

s 

1 

* 1O -0 

5 

He2 

♦ 

2e ■* He* ♦ He ♦ e 

t 

(4 t 0.5) 

x 10’ 2 ° 

6 

HeJ 

♦ 

e ♦ He He* ♦ 2He 

u 

(5 ± 1) 

x 10- 27 

6 

He* 

♦ 

He* -► He 2 * e 

X 

(1.5 ± 0.3) 

* 10- 9 

6 

He* 

♦ 

e ♦ He ♦ e 

y 

(4.2 ± 0.6) 

x 10- 9 

6 




Table 4.2. 

Argon rate 

coefficients. 




Reaction 

Symbol 

Rate Coefficient 

Reference 

Ar* 

♦ 

2e -*■ Ar* ♦ e 

t' 

~ 9 x 10 -20 

7 

Ar* 

♦ 

e ♦ Ar ■* Ar* ♦ Ar 

u' 

~ 1 x io -26 

7 

Ar* 

♦ 

2Ar ■* Ar 2 * Ar 

w ( 

>•*:!:?)■ io - 31 

8 

Ar2 

♦ 

e -► Ar* ♦ Ar 

*2 

(7 ± 1) x 10- 7 

9 


Table 4.3. Helium-argon rate coefficients. 

Reaction Symbol Rate Coefficient Reference 

He 2 ♦ e ♦ Ar ■+ He* ♦ He ♦ Ar ~ u 5 * 10"*- 7 

He* ♦ Ar -*• Ar* ♦ He ♦ e p (7 ± 1) * 10” 11 10 

He£ ♦ Ar Ar* ♦ 2He q 2.2 * 10" 10 11 

He 2 ♦ Ar ♦ He ■* Ar* ♦ 3He r (24 ± 6) * 10" 30 11 

Ar* ♦ Ar ♦ He •♦ Ar 2 ♦He v ^ 1 * 10" 31 8 

Ar* ♦ e ♦ He - Ar* ♦ He ~u' -1 * 10" 26 
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Figure 4.1. Flow chart of processes in 
the 3 He-Ar Laser. (Dashed 
lines show source terms). 
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Figure 4.2. Comparison of approximate analytical solutions for the 
densities (dashed) and exact computer solutions 
(continuous lines). 
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Figure 4.3. Processes occurring after recombination. 

The processes shown with dashed lines are 
unimportant . 
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Figure 4.4. Comparison with experiment. 
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